Exchange effects in plasmas: the case of low-frequency dynamics 
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Recently, there has been a surge in the interest of non-equilibrium collective quantum models, where particle 
dispersion and spin are examples of effects taken into account. Here, we derive a kinetic plasma model con- 
taining fermion exchange effects. Exchange interactions are of great importance in many systems, and have no 
classical analogy. Our model therefore constitute a possible probe of collective quantum phenomena in new 
regimes. As an example, we consider the influence of exchange effect on low frequency dynamics, in partic- 
ular ion acoustic waves. Comparisons to related computational techniques are given and the differences are 
highlighted. Furthermore, we discuss the applicability of our model, its limitations and possible extensions. 

PACS numbers: 52.25.-b, 52.25.Dg 



43 ' 

or 

I 



ll- 
CC '. 

O ■ 

• i-H , 

C« . 

>v 

Of 



> 

ON 
00 

m 
cn 

cn 
o 



x 



Quantum plasma physics is currently a field of intense 
study. One reason for this is the potential applications in, 
for example, laser produced plasmas, ultra small electronic 
devices, and dense astrophysical systems (Urfl. Different as- 
pects of quantum plasmas have been studied such as quan- 



3], the magnetic dipole 
quantum relativistic ef- 



1J, LLJ- Typically, 



turn dispersion and Fermi pressure 
force and the spin dynamics 14413 
fects and nonlinear dynamics fllfj|-ll2 
quantum effects are important for systems with high density 
and low temperature. This said, it is important to distinguish 
between quantum effects related to thermodynamic equilib- 
rium properties and dynamical properties of the system. Ex- 
change effects due to particle statistics have been successfully 
included in the density functional theory (DFT) 11161 Il7ll . Ap- 
plications of DFT include for example ground state properties 
of atoms and equilibrium properties of many-particle systems 
via density functional theory Ill6lll7ll . The effects of exchange 
on dynamics have also been studied in the setting of kinetic 
theory ii, 18-24], as well as in studying, e.g., the thermody- 
namic properties of plasmas 1 2(1 21 1. It has also been studied 
using fluid theory Il24ll . Furthermore, many papers deal with 
how quantum mechanics affects the low-frequency long-scale 
dynamics, as for example quantum ion acoustic waves ]25 



In the present letter we have derived the quantum kinetic 
evolution equation for electrons within the Hartree-Fock ap- 
proximation, treating the exchange effects perturbatively. The 
general result is further simplified by considering a plasma 
that is not spin-polarized, and by focusing on length scales 
much longer than the characteristic de Broglie length. We ap- 
ply the formalism to one-dimensional ion-acoustic waves in 
plasmas, within the linear approximation. 

We here consider a completely ionized electron-ion plasma 
with the particles interacting through a mean-field scalar po- 
tential. Quantum effects for the ions will be completely ne- 
glected, while for the electrons we will take into account a 
dynamic correction due to the Pauli exclusion principle. Also, 



we will not consider effects due to the self-energy and particle 
correlations ll2l1l . We will now give an outline of the deriva- 
tion of a kinetic theory, where the exchange effect is taken into 
account. 

The state of the A^-electrons is described by the density op- 
erator pi. jv (see for example Ref. |2lf|). and the dynamics is 
given by the von Neumann equation with the Hamiltonian 
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Here m e is the electron mass, e is the electron charge (<? < 0) 
and £o is the permittivity of vacuum. The last term accounts 
for the interaction with the electric potential created by the 
ions. We now introduce the reduced density operators ac- 
cording to pi i = «'Tr, + i jvPl.jvAi...;, where n is the mean 
density and Ai , is the operator that takes an /-particle state 
and makes it completely antisymmetric. We will only need to 
know that A12 = 1 — P\i where Pyi interchanges particle 1 and 
2 (see, e.g., Ref. J2lll for further details). The evolution for 
the one-particle density operator is given by 



ihdtPi = [/!i,Pi]+«Tr 2 [Vi2,pi2Ai 2 ], 



(2) 



where h\ = p / (2m e ) and V12 = V{k\ — £2) = q I (4ft£o|xi — 
X2I) and P12 is the two-particle density operator. The effects 
of two-particle correlations gn can be separated out of the 
two-particle density operator by writing it in the form P12 = 
P1P2 +gi2 ■ We are interested in the collisionless limit where a 
mean-field approximation will suffice. This approximation is 
obtained by neglecting the correlation function gi 2 . Utilizing 
this in the equation above we then get 



;Mpi = [/ii,pi] + [Vi,Pi] 



(3) 



where V\ = Tr2Vi2P2Ai2i is the Hartree-Fock potential. This 
is a closed system for the one-particle density operator. The 
evolution equation Eq. ([3]l can also be written in the Wigner 
representation, i.e. 
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where 
/(x,p,a,/3) 
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ternal) potential, respectively and V(x) = ^ 2 /(47T£o|x|) is the 
Coulomb potential. The LHS of Eq. represents the quan- 
tum Vlasov equation, while the RHS is the correction due to 
exchange effects. This term is nonlocal in phase-space and 
(5) nonlinear in the distribution function. 



is the Wigner matrix for a spin-1/2 particle (a,j3 = 1,2, for 
spin up and down along the axis of quantization) and 



The matrix equation can be transformed into a scalar equa- 
tion by taking the spin transformation JH 



0(x) 



/(x,p,s,/) = — £ [S aiP +s-a atP ]f(x,p,P,a), (7) 



a,B=l 



where s is a vector on the unit sphere. Applying this to Eq. 
is the Wigner distribution and the total (mean-field and ex- we obtain 
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where in the last term we see the exchange interaction in the 
Wigner form. 

The evolution equation describes the evolution of the 
electrons in the mean-field approximation for all scale lengths. 
We are interested in the semiclassical limit where the poten- 
tial and the distribution function / vary on a scale L much 
larger than the de Broglie scale length A^b and would like to 
only keep the lowest surviving correction in an expansion in 
Aaq/L. For the potential term the expansion is straightfor- 



l x + r , y „' A f \ x + r y „" c" 



(8) 
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ward, see for example Ref. For the exchange term we 
expand the potential and the distribution function to second 
order in y (with the assumptions that the characteristic length 
scale L is much larger than the thermal de Broglie wave length 
h/rnvj, where vj is the thermal speed). We then perform the 
y-integration and one of the momentum integrals. 

Furthermore, we will for simplicity also assume that 
the distribution function is independent of the spin, i.e. 
/(x,p,s,f) = f(x,p,t)/(4n). Integrating over the spin we ob- 
tain 
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where d l x = d/d x, and analogously for d' and an arrow above an operator shows which direction it acts. We have also 
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used the summation convention, so that repeated indices are 
summed over. 

We now consider the effect of the exchange term on elec- 
trostatic ion-acoustic waves in a plasma. We will use Eq. 
© for the electrons and the classical Vlasov equation for 
ions. To obtain the dispersion relation we assume a longitu- 
dinal oscillation / = fo(p) + f\ (p)exp(— icot + ikz) and E = 
z£ z exp(— icot + ikz). We assume that the unperturbed elec- 
tron distribution function is given by a Maxwell-Boltzmann 
distribution H 

/Q(P) ^ (2^)3/ 2 eXP ('^l)- (10) 

Furthermore, assuming that the exchange terms are small cor- 
rection to the distribution function, we may calculate it by in- 
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serting the lowest order solution for f\ in the integrand. In- 
troducing spatial spherical coordinates, it is possible to solve 
all spatial integrals if one assumes that the integrand vanishes 
in the limit r — > °°. Next, the integrand is expanded in terms 
of h. The lowest order term in the first integral vanishes due 
to symmetry and we keep only the two first order terms. In 
the second integral we already have an additional h, meaning 
that we only retain the lowest order term. Finally it is possible 
to solve the pi and p' integrals. The remaining integrals are 
solved numerically and doing so gives a solution for f\ in the 
linear regime. Now, from the classical dispersion relation we 
have 

co « {(O pI /(O pe )kv T e = akv Te , (11) 

where vr e = ^/kBT e /m e is the electron thermal velocity. The 
dispersion relation is then given by 
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where the first three terms gives the classical dispersion rela- 
tion for an ion-acoustic wave. Solving these integrals numer- 
ically gives the approximate dispersion relation 



0« 1 



co- 



k 2 v\ e 



2/y c i 

kc s 



co 



Pi 



co 2 m 2 k 2 \'j e 

2,2 



(0.8+0.05/). 



which in the quasi-neutral limit ©£, ;§> k l vf e can be written 



co = kc s \ 1+0. 
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where c s = (m e /mi) l ' 2 VTe is the classical ion-acoustic ve- 
locity and we have introduced the classical electron Landau 
damping, jci = kc s ^7t/8^m e / m,-, in the cold ion limit 131 1. 
Thus we note that the effective ion-acoustic velocity is in- 
creased, whereas the damping due to wave-particle interaction 
is decreased due to the exchange effect. As seen from ([TBI 
the relative magnitude of both these effects is of the order H 2 , 
where H = hco pe /m\'j e . As is shown in, e.g., Ref. HI, plotting 
the line H = 1 in a log-log density temperature diagram di- 
vides the parameter space in a classical regime (H 1) and a 
strong quantum regime (H > 1). However, such plots are typ- 
ically performed in order to illustrate the relative importance 
of particle dispersive effects. Within a kinetic formalism par- 
ticle dispersion is described by the terms with higher order 
momentum and spatial derivatives in the Wigner equation 111 . 
For such terms to be important, in addition to the parame- 
ter H not being too small we also require the macroscopic 
scale lengths under study to be short. Specifically we need 
the scale lengths to approach the thermal de Broglie wave- 
length or shorter. Thus if we exclude the regime of short scale 



lengths, as we have done here, the quantum effect of particle 
dispersion is guaranteed to be of little significance. By con- 
trast, we see that exchange effects may very well affect the 
long scale behavior of the low-frequency density dynamics. 
Of particular interest is the change in the damping term. By 
approaching the regime H ~ 1, Eq. dot suggests that we may 
more or less completely suppress Landau damping of ion- 
acoustic waves. Physically this makes sense, as classically 
the particles that are resonantly accelerated for a long time 
are rather well localized in phase space, which is then coun- 
teracted by the exchange terms. Strictly speaking, the regime 
H ~ 1 does not fit into the perturbation scheme that we have 
applied here, but qualitatively we still expect this result to be 
valid. 

An important result from this study is the general expres- 
sion for the exchange term, as given by Eq. (|9}- This term can 
describe exchange modification of any type of processes, e.g . 
altering the coefficients for three-wave interaction 11321 13311 . 
adjusting the Zakharov equations D4 13511 or modifying non- 
linear wave particle interaction processes 113 6[ |37fl . The main 
restriction is due to the assumption of electrostatic fields. The 
complexity of the exchange interaction term in ((9) in practice 
forces one to do perturbative calculations. Since the present 
formalism captures the full effect of a distribution function 
which may be far from equilibrium, it provides a valuable op- 
portunity to evaluate approaches that relies on other types of 
approximations. Specifically, in time-dependent density func- 
tional theory (TDDFT) 0811 the properties of the system is 
derived from the electron density only (or is at least limited to 
macroscopic quantities), in which case the full dependence on 
the detailed momentum distribution is disregarded. Due to the 
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complexity of many nonlinear plasma systems, such a drastic 
simplification may be needed, but at the same time is it essen- 
tial that the accuracy of the approach can be evaluated. Re- 
sults from DFT calculations have been used to describe elec- 
trostatic waves in plasmas, see e.g. Eq. (6) of Ref. J39ll . 
where the further approximation of the adiabatic local density 
approximation (ALDA) has been used. However, a difference 
with our case is that the Fermi temperature was assumed to 
be higher than the plasma temperature in these papers. In a 
very rough sense the previous results agree with ours, as the 
relative importance of the exchange term scale as (Tkq p /Ek) 2 
in both cases, noting that the characteristic kinetic energy 
is the thermal energy Ic^T in our case and the Fermi energy 
kgTp in the case of Ref. 13911 . However, in our case the phase 
velocity of the ion-acoustic waves is increased due to the ex- 
change interaction, whereas based on Eq. (6) of Ref. 0911 the 
phase velocity is decreased. Still the interpretation of this fact 
can be debated. One possibility is that the approximation of 
ALFA to evaluate the exchange potential is too restrictive to 
capture the ion-acoustic dynamics accurately. Another possi- 
bility is that the results are indeed sensitive to the ordering of 
T and 7>, such that the sign of the exchange effect is reversed 
when the ordering is changed. Regardless of this, is is clear 
that DFT calculations in general cannot capture the effects of 
wave-particle interaction , which is responsible for the wave 
damping in our case. 
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